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I. INTRODUCTION 

Evecy day, military aad civilian forecasters are tasked 
to provide weather forecasts for remote areas of the world. 
Many are short-range forecasts (up to 12 h) for specific 
locations, rather than large areas. Facsimile charts 
combined with climatology provide reasonable wind and temp- 
erature forecasts, which can be cross-checked using satel- 
lite imagery. in addatvon, satellite imagery aids in 
preparing sky coverage forecasts. A major remaining problen 
is precipitation. Since there may be few, if any, timely 
surface observations available, 1t can be Zgafficult Was 
determine whether precipitation is falling beneath areas of 
cloud cover. What is required are metnods to differentiate 
precipitation from non-precipitation uSing geostationary 
visual and infrared satellite data. 

Research in this area iS important for two inter- 
changeabie reasons. First, better okservations will lead to 
better forecasts. A fundamental problem with all fore- 
caSting schemes 1s a lack of reliable initial data. An 
accurate and reliabie scheme for differentiating between 
precipitating and non-precipitating clouds will greatly 
improve current and future forecasting schemes. Second, the 
use of satellite data will increase the area of coverage and 
therefore increase the number of potential observations. 
The more "good" initial data the better the chance the fore- 
castiny scheme will have to produce a better forecast. 

Although Microwave sensing has shown much promise int 
detecting rainfall, visible and infrared data are utilized 
in thie reseanems A major drawback of microwave sensing 1S 
that it uses a broad field of view and has a more coarse 
resolution. In addition, visible and infrared data are 


currently more readily available (every 30 minutes). 
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Previously, Studies have been undertaken which have 
attempted to correlate precipitation with satellite imagery. 
Muench and Keegan (1979) studied GOES-E visual and infrared 
Satellite imagery and hourly rainfali climatological data 


momemrave surface Stations in the northeastern United States 


for the period April through November 1977. Their data set 
consisted of 2760 surface observations, of which 300 
reported precipitation. Using infrared and visual imagery, 


they attempted to predict hourly rainfall amounts. Although 
their results were "far from ideai", they were still "useful 
for making decisions if the threshold for ‘action' versus 
‘no action was known" (Muench and Keegan, 1979). 

In addition to Muench and Keegan i197 Fie other 
researchers have studied this problen. ibaisas- (196 lar, 
1981b) developed a bi-spectral cloud classification based on 
@iemar and infrared data from the polar-orbiting TIROS-6 
satellite. Loves oy -ancwhiictinm, (19798, using the results of 
Muench and Keegan (1979), developed a two dimensional 
freyuency piot for cumulus and non-cumuluS precipitation. 
In addition, they also tested a Spectral threshold techniyue 
for rain area mapping. For those unfamiliar with research 
in this area, a summary by Paul (1983) of selected studies 
is included in Appendix A. 

mae iges), using only daylight observations (0800L - 
1600L EDT), computed means ard standard deviations for each 
mo x 10, Se oP nOam Sora naed. “el Matrix Of pixels (both 
infrared and visual) centered on a surface reporting station 
with collocated geostationary satellite data. Her results 
Showed a distinct separation between the means cf the 
precipitation and non-precipitation cases. However, adding 
one standard deviation to the infrared precipitation mean 
and subtracting one Standard deviation fron the non- 
precipitation mean resulted in a Significant overlap between 


precipitation and HONM=precipitation thresholds. os 
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overlap was too iarge to permit an unbiased forecast to be 
made. Similar results were found using visual sateliite 
data. 

This research will concentrate ona relating sateliite 
imagery to observed precipitation and fnon-precipitation 
events using the surface-satellite data set of Pau1l(1983). 
The objectives are: 

1. Use a full 24 h data base, vice the 9 h data Lase by 
Pawel ie3) 3 

2. Form distributions of satellite data for several 
classes of observed weather events, test for normality ane 
study the separation of tke distributions. 

3. Apply established visual and infrared thresholds to 
the data set and determine their reliability. 

This thesis is organized into four chapters. Chapter II 
describes the data set. Chapter III presents an analysis 
and discussion of the data _ set. Selection of matrix sizes 
for analyses and choices for studying the characteristics of 
each distribution are explained. Results are then compared 
to Paul (1983). The £ ide section verifies the 
precipitation/no-precipitation infrared and visible thresh- 
olds currently used in SPADS. Chapter IV states conclusions 


and suggests further study. 
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if. DATA DESCRIPTION 

The data set assembled for this study consists of collo- 
cated GOES-E satellite data and Service-A hourly surface 
observations for the southeastern United States during 
magust 1977. During August, this region is dominated by 
Subtropical airmasses with extensive convective activity. 
ihe, GOES-E.data consists of 10 x10 »bixel matrices of 
visible and infrared satellite data centered over each of 
137 reporting stations (Fig. 1), all south of 4O N. The 
Satellite data are measured with the Visual Infrared Spin 
Scanned Radiometer ia SSR) whick have subsatellite point 
Spatial resolutions of 1 km for visuai channels and 7 km foc 
infrared channels (U. S. Department of Commerce, 1983) The 
GOES-E navigation was completed by Man-computer Interactive 
Data Access System (McIDAS) at the University of Wisconsin 
feamg the full resolution visual data, with an accuracy of 
lemepixels (1-2 km). The full resolution visual data were 
averaged to a 7 km resolution, to egual the infrared data 
resolution. The visual and infrared digital counts ranged 
mean Values of 0-63 and 0-255, respectively. The 10 x 10 
GOES-E visual and infrared satellite data each cover an area 
Cron mi x 45 nmi at SO N (60 nmi x 60 n mi at 42 N). 

What distinguishes this data set from those used in most 
Other studies is the number of obServations. This data set 
Naoea total of 70,623 Service-A hourlyereports, o£ which 
29,342 reports have collocated visual and/or infrared satel- 
lite imagery. The visual data were normalized and converted 
to albedos based on the algorithm of Muench ana Keegan 
(1979). Their scheme corrects for the varying zenith angles 
as well as adjusting the visual satellite data for aniso- 


tropic scattering as related to the zenith angle. While 
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albedo values cannot exceed 1.00, this scheme allows vaiues 
to overshoot 1.00 up to a value of 1.20. Therefore, the 
visual satellite values are not true albedos, but are esti- 
mated albedos. The extended visual normalized data scale 
was used to permit comparison of the results of this 
research to the bi-spectral threshold specifications of 
Muench and Keegan (1979) and to Paul (1983). The Muench and 
Keegan normalization scheme specifies that any computed 
albedo greater than 1.20 be set egual to 171.20 to Jlinit 
unreasonably large values. Similarly, the scheme s;,ecifies 
computed albedos less than 0.15 be interpreted as the ground 
or water surface reflectance and the value 0.00 be assigned. 
(See Appendix B for further specific information concerning 
the Muench and Keegan normalization scheme.) The infrared 
data were processed in digital counts and converted to cloud 
top temperatures prior to statistical calculations. 

The surface reports within the data set contain informa- 
tion that regularly appears in ail Service-A reports, plus 
additional location information and viSual and iNfrapged 
satellite data (see Table I). Of the information available, 
we were most concerned with the current weather, the cloud 
group, and the satellite data. 

For this study, the data set was divided into four 
categories: 

1. Observations with ho precipitation and inftfaned 
Satellite data (19,354 Obs) ; 

2. Observations with no precipitation and visual satel- 
lite data (14,740 Obs); 

3. Observations with precipitation and infrared satel- 
lite data (1856 Obs); and 

4, Observations with precipitation and visual satellite 
data (1437 Obs). 

The no-precipitation data category iS comprised of all 


observations not showing any "R" ain tne current weather 
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mou p. item beaGcittdtaon Gdatd Category consists or all 
Seeervations with an “R" in the current weather group. 

The cloud cover group 1S based ona three digit code. 
The first digit indicates the amount of low clouds, where 0 
is defined as clear, 11s scattered, 2 is broken, and 3 is 
overcast. The definitions for clear, scattered, broken, and 
overcast are as defined in the Federal Meteorological 
Handbock No. 1 ({U. S. Department of Commerce, 1980). The 
second and third digits of the cloud group indicate the 
amount of middle and high clouds respectively. 

In this study there were no restrictions to a specific 
time frame. Rather, all observation, regardless of time, 
which had satellite data available were used. The only 
restrictions imposed were: 

1. Only broken or overcast skies were considered; and 

2. Each pixel count must be greater than zero. 

Only broken or overcast ceilings were chosen to be analyzed, 
based on the assumption that precipitation from scattered 
skies is rare and that precipitation from clear skies iS an 
MacoLrrect observation. indeed, Paui (1983) reported that 
less than one percent of the precipitation observations 
within her time period also reported scattered skies. The 
restriction that each pixel count be jgreater than Zero was 
imposed because any fixel less thanh zero in meauingless. 

Paul (1983) limited her time period to daylight hours 
many (0800 —- 1600L EDT) to avoid distortio: of the visible 
data due to low solar elevation angles. By choosing not to 
jimit the time period for this study, it was understood that 


HrSslLOrtL1Oon Of visible data may occuL. 
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III. DATA ANALYSIS 


A. INTRODUCTION 


This chapter consists of analyses and discussions of the 
data set. It begins with a discussion on choosing the most 
representative and reliable matrix size and the distribu- 
tions formed by those matrices. Next, the distributions are 
tested for normality. Using those results, the distribution 
characteristics are described and discussed. Next the 
results are compared to those obtained by Paul (1983). In 
the final section, the data set is used to verify the estab- 
lished precipitation/no-precipitation thresholds of an auto- 
Mated cloud and precipitation estimation scheme (Wash et al, 
1984) with the U. 5S. Navy's Satellite Data Processing and 
Display System (SPADS). 


Be MATRIX SIZE 


The collocated satellite data consisted of a 10 x 10 
matrix centered on the surface observation station. This 
provided the option of using any one of nine different 
Matrix sizes, 10 x 10 through 2 x 2, in the research. Paul 
(1933) evaluated four sizes, 10 x 10, 8 x 8, 6 x 6 and 4 x 4Y 
(2925 sq nmi, 1296 sG nmi, 729 Sg nh Bi and 46% Sq Hime 
respectively). Her results (using 7358 no-precipitation 
observations and 534 precipitation observations) showed no 
Significant differences in the neans and standard deviations 
of the four matrices. When Paul's time restrictions were 
removed, the results (using up to 19,354 no-precipitation 
observations and 1856 precipitation observations) also 
showed no significant differences between the same four 


matrices based on the means and standard deviations (Table 
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Zi). From these results it was concluded that the distribu- 
tion of the data would not change significantly as the 
matrix sizes were changed. The matrix size could now be 
based on meteorological parameters. 

Initially one could argue that the 2 x 2 matrix would be 
the most representative. It could most accurately represent 
the surface observation by not being influenced as much by 
the surrounding area as any of the larger sizes. The major 
disadvantage of the 2x 2 matrix though, is actually its 
Smail size. Because the navigational error of McIDAS is one 
to two pixels, the possibility that the surface observation 
would actually be outside tke satellite data matrix is 
present. Choosing the 4 x 4 matrix would greatly reduce the 
possibility of having the surface observation outside the 
Sateilite data matrix. The 6 x 6 matrix would guarantee 
fWathin McIDAS error) the surface observation being within 
the matrix, but now areal coverage must be considered. The 
Gex 6 matrix has an areal coverage of 729 sq n mi at 30 N 
while the 4 x 4 matrix covers 484 sq n Mi, approximately 66% 
@™ the 6 x 6 matrix. After weighing the possibility of a 
misplaced observation using the 4 x 4 matrix, with the 
increased coverage of the 6 x 6 matrix, the 4 x 4 matrix was 


chosen as the most representative matrix siZe. 


C. NORMALITY 


The determination of the hormality of the distributions 
within the data set waS an important stipulation of this 
research. Since the normality determined the way in which 
the distributions could be described (using confidence 
devels if normal, and just general descriptions if not 
normal), it waS paramount that the distribution be tested 


SOrLect ly. 


ed 


Histograms were plotted for each of the four data 
categories with observations which had collocated satellite 
data. Within each category; ea 10° x90, 8 x 38, Gia) andi 
4% matrix was also plotted. Even though the 4 x 4 matrix had 
already been chosen as the size to be evaluated, the option 
to change sizes if warranted by the normaiity testing was 
left open. 

Toe Chi-square test was chosen to test the distribution 
fOr neLialaty. It involves comparing a sample distribution 
to a normal distmapucion. Each distribution was divided 
Into.) 20 Seba: Using the mean and standard deviation 
computed earlier for each category, 17 degrees of freedom 
were available. All of the no-precipitation categories 
could be seen not to be normally distributed by inspection 
of the nistograms (see Fig. 2 and Fig. 3) and therefore were 
not formally tested. Since each precipitation category nore 


closely resembled a normal distribution, they were tested 


with the Chi-scuare test. Each no-precipitation distribu- 
tion tested (all four matrix sizes) proved not to be 
distributed normally. Fig. 4, although appearing to be 


normally distributed, only had a 5% chance of coming from a 
normal distribution, and it was the closest to normai of any 
category tested. 

It must be understood that just Decause a distribution 
kas a low confidence level for being normally distributed, 
it does not mean that it is not anormal distribution. 
Likewise, just because a distribution has a high confidence 
level of being normally distributed, does not mean that it 
is a normal distribution. Since the results Showed a nigh 
probability of not being normally distriktuted, the research 
focus shifted to describiny the non-normal and other inter- 


esting aspects of the distributions. 
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fee DISTRIBUTION CHARACTERISTICS 


- Int roduwetrion 


This section will present a discussion about each of 
the histograms plotted for infrared/visual and 
no-~precipitation/precipitation. It is important to remenber 
that only broken or overcast ceilings are included in the 
Sestributions. LisOpadulcu™eomLrouna-off™ error, not all 
percentages calculated in this section add up to 100% 


exactly. 


2. No 


This category consists of all surface observations 
without precipitation and with collocated infrared sateiiite 
data (see Fig. Die There were 19,354 observations within 
this category. As expected, the majority of the observa- 
tions are grouped near the warm end. .Tae cioud top tempera- 
tures ranged from a @Binimum of 205.0 K to a maximum of 320.1 
Ne It has a peak of 4214 obServations between 290.0 K and 
er 0 K. The mean temperature was 275.7 K with a standard 
deviation of 19.4 K. 

The most interesting part of this distribution is 
the cold tail (less than 240.0 K). These are observations 
which would ke misclassified by a single cold threshold. It 
was found that of the 1461 observations within the cold 
tail, overcast and broken ceilinys were divided evenly (see 
meavode III). Within the overcast ceiling category, 35% were 
high clouds with 10% mid-level clouds and 5% low clouds. 
Within the broken ceiling category, 41% were high clouds 
with 7% mid-level clouds and 3% low clouds. If two or more 
levels reported broken ceilings, tne highest ceiling was 
specified as the ceiling for this study. For example, if a 
station reported 3500 broken, 10, 000MbrEOken; and 22, 000 


broken, the ceiling would be reported as a high ceiling. 


US 


Combining the overcast ceilings with the broken ceilings 
produced 76% of the okfservations with high clouds, 16% with 
mid-level clouds, and 8% with low clouds. 

The breakdown of the overcast ceiling category can 
be misleading. When an observer reportS an overcast 
ceiling, he cannot see what is above the base of the clouds. 
The misrepresentation arises because the tops of the clouds 
are not at the level reported py the observer. In addition, 
there may be additional layers above the overcast layer. 
The satellite data wili therefore contain information about 
the highest (coldest) Clouds it senses and that layer may 
possibly he one or two layers above the reported ceiling. 

TO approximate the percentage of overcast ceiiings 
affected, the high cloud with broken ceilinys category was 
further analyzed. It was found that 95% of the observations 
within that category were multi-layered clouds (with either 
scattered or broken conditions reported at low or nid-levels 
in addition to a broken ceiling reported at the high level). 

The overwhelming nuwber or warm cloud top tempera- 
tures 1S not surprising for a non-precipitation category. 
Of some concern though, are the few temperatures (66) which 
exceed 300.0 K (27.0 C). One explanation is that the 
infrared sensor is seeing the top of a marine layer along 
the Gulit coast. Another possibility is that the surface of 
the earth is influencing the mean vaiue of the 4 x 4 matrix. 
This is possible because one surface observation is repre- 
senting a 424 sq n miarea. Just because a ceiling is being 
reported at the observation location, does not necessarily 
mean there 1S a ceiling 19 n mi away. Warm surface tempera- 
tures could increase the mean value of the matrix 
Signitieantay. 

in summary, the large majority of warm temperatures 
for non-precipitation cloud tops 1s just as one would 


expect. The extremely warm temperatures do not affect the 
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Mecipitation/no-precipitation decision. iieswco1.d tail is 
composed mostly of thick high clouds or overcast/broken 
Clouds with additional clouds at lower levels. Signe icant 
errors arise if these cold cloud tops are interpreted as 


Mgeecipfpitating clouds. 


This category consists of all observations without 
precipitation and with collocated visual satellite data (See 
Fag. Bs There were 14,740 observations within this 
category. As expected, the majority of the observations are 
grouped near the dim end. The estimated albedos ranged from 
meee, to 1.20. It has a peak of 1944 observations at 0.00. 
The mean albedo was 0.27 with a standard deviation of 0.26. 

Further analyses were conducted on the 999 observa- 
meons COMprising the bright tail (greater than 0.70). which 
would lead to misclassification of precipitation from the 
wecwal data. It was found that 49% of the observations 
reported overcast ceilings and 51% reported broken ceilings 
(see Table IV). Within the overcast ceiling category, 18% 
were high clouds, 14% were mid-ievel clouds, and 17% were 
iow clouds. Within the broken category, 41% were high 
clouds, 7% were mid-level clouds, and 3% wece low clouds. 
Combining the overcast ceilings with the broken ceilings 


produced 59% of the observations with high clouds, 22% with 


mid-level clouds, and 20% with low clouds. As in Section 
[eeeZ , the breakdown of the overcast ceilings can be 
misleading. TO approximate the actual percentage of high 


clouds, further analyses were conducted on the high cloud 
with broken ceiling category for multi-layered clouds. et 
was found that 83% of those observations contained multi- 
layered clouds. 

The large number of albedos between 1.15 and 1.29 is 


@ause for concern. Paul #1983) , by imposing a time 
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restriction on the data base. excluded most of the problems 
associated with distortion caused by low solar angles near 
sunrise and sunset. Because this study did not restrict the 
data base to a specific time period, distortion waS a 
problen. In the distribution for this category there were 
336 reports with albedos greater than 1.15. Further anali- 
ysis of these reports showed all albedos greater than 1.15 
were reported at 0000 GMT (near sunset). In addition, ar, 
of all observations with albedos greater than 0.70 occurred 
within three hours of sunrise or sunset. In addition 
high albedo values caused by distortion or nearby precipita- 
tion, further analysis Showed that one-third (44/135) of the 
observations betweem 0.70 and 1.15 which reported low over- 
cast ceilings also reported fog. It must be noted that 
although these observations were clustered around sunrise, 
they were still assigned higher albedo values than one would 
normally expect for fog. A significant portion of these 
bright non-precipitation observations are related to prob- 
lems in normalizing the satellite data for low sun angle 
Situations. More work iS reguired in using the cloud 
brightness normalization before more quantitative use of the 
data can be made. 

In summary, most of the bright tail is composed of 
thick high clouds or overcast/broken clouds with clouds also 
at lower levels. AS in the no-precipitation, infrared 
category, Significant errors arise when these bright clouds 
are interpreted as precipitating clouds. In addition 
distortion due to low solar angles also causes misclassified 


bright albedo values. 
4. Precipitation, Infrared 


This category consists of all observations reporting 
precipitation and with collocated infrared sateilite data 


(see Fig. 5). °(There were 1856 observations within this 


a 


Gategory. The distribution of cloud top temperatures (clus- 
tered near the cold. end) 1s Similar to what would be 
expected for a precipitation category. The cloud top temp- 
eratures ranged frora minimum of 205.0 K to a maximum of 
moioe 6 K. It has a peak of 181 observations between 220.0 K 
ame 225.0 K. The mean is 243.0 K with a standard deviation 
@i 21.7 K. 

Further anaiyses were conducted on the warm tail 
(warmer than 270.0 K), the cold end (colder than 240.0 4X), 
and the area between them {see Table V). The warm tail was 
analyzed because it would most likely contain misclassified 
observations. The cold end and the middle area each were 
analyzed to better understand the satellite data distribu- 
ON . Each of the three areas were analyzed for the type 
and intensity or the precipitation, and zor convective 
versus stratiform precipitation. 

In the warm tail there were 257 observations. Light 
rain showers were dcminate with 50% of the observations. 
megane rain and light thunderstorms followed in frrequency 
with 21% and 18%, respectively. An overwhelming majority, 
91%, reported light precipitation, compared to only 9% with 
moderate/heavy precipitation. in. Gaddit ion; 76% of the 
observations were convective in nature compared to 24% of a 
stratiform hature. One can expect the observations within 
the warm tail to be misclassified becaise precipitation is 
not usually associated with warm cloud top temperatures. It 
Should be noted though, that only 9% of the precipitation 
observations within the warm tail reported moderate or heavy 
PEeGipitation, whereas 91% reported light precipitation. 
This leads one to conclude that although the observations 
within the warm tail could be misclassified, they would 
would be light precipitation and not moderate/neavy 


observations. 


23 


A possible reason why some pireClipatatioe 
precipitation observations have Such warm cloud =tops me 
because the precipitation, in addition to being light, is 
also very scattered. In this situation, ©a »Oort1on ofee 
satellite data is composed of surface reports rather than 
cloud "cor edata. Also, it should be noted that the surface 
observations are taken 5-10 minutes before the hour, while 
the satellite Scan is taken 5 minutes after the hour. Thais 
10 to 15 minute lag in the satellite scan 1s at least 


partially responsible for some of the misclassified observa- 


tion's < Precipitaticn could be reported at the surface 
station at the observation time, but by the satellite scan 
time, it may have stopped. Tnais would produce a surface 


precipitation report with a collocated ho-precipitatiag 
satellite report. 

In the cold end there were 959 observations. Ligne 
thunderstorms were reported 42% of the time with light rain 
Showers being reported 20% of the time. In addition, 15% of 
the observations reported moderate/neavy thunderstorms and 
14% reported light rain. There were more light precipita- 
tion reports (81%) than moderate/neavy reports (194). Once 
again convective type precipitation (79%) was reported more 
often tnan stratiform precepitation (Zig 

In the region between the warm tail and the cold end 
there were 640 observations. Light rain showers were in 37% 
of the observations, followed by light rain in 29% and light 
thunderstorms in 21% of the observations. Light preci pita— 
tion was reported in 89% of the observations, with moderate/ 
heavy precipitation in the remaining 11%. As before, 
convective type precipitation (66%) waS reported more often 
than stratiform precipitation (34%). 

It should be noted here that this study cannot 
distinguish between precipitation intensities with any 


degree of accuracy. One would expect a much higher 
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percentage of moderate/neavy thunderstorms in the cold end 
than in the warm tail. These results snowed only a slight 
increase (9% to 19%) in occurrence of heavier precipitation 
moving from the warm tail to the cold end. The aajor reason 
is that this study uses an average of surrounding points to 
generate a collocated satellite observation, whereas other 
studies, which show a distinct increase in precipitation 
intensity with colder cloud top temperatures, used the 
coldest cloud top temperature as the value. 

In summary, although the occurrences oz moderate/ 
heavy precipitation were less than expected, tne trend of 
increaSing intensity with colder cloud top temperatures was 
present. Aliso, the trend of increaSing thunderstorms and 
decreaSing rain showers with colder cloud top temperatures 


is consistent with expectations. 


5. Precipitation, Visible 


This category consists of all observations which 
reported precipitation and had collocated visible satellite 
data (see Fig. 6). There were 1437 observations within this 
category. The distribution 1S aS one would expect fora 
precipitation with visibie satellite data category (albedo 
values clustered near the bright’ end). The estimated 
albedos ranged from 0.00 to 1.20. There is a peak of 132 
observations between 0.65 and 0.70. The mean albedo is 0.61 
Math a standard deviation of 0629. 

This distribution was divided into three areas and 
further analyses were conducted on the dim tail (albedos 
less than 0.40), the bright tail (albedos greater than 
0.90), and the region between the dim and the bright tails 
(see Table VI). The dim taii was chosen because it 1s the 
area of the greatest probability of misclassified observa- 
tions. The other two areas were chosen to compare to the 


dam tail. Within each of the three areas, an analysis for 
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precipitation type, precipitation intensity and convection 
or stratiform precipitaervon was conducteur 

In the dim tail there were 327 observations. Li Gate 
rain showers (43%) were dominate with light thunderstorms 
(28% and light rain (15%) following in freguency. Light 
precipitation was reforted in 88% of the observations with 
moderate/heavy preci;fitation in 12% of the observations. [In 
addition, convective type precipitation was reported 81% of 
the time with stratifcrm precipitation being reported 19% of 
the time. As expected, most of the misclassified observa- 
tions were light precipitation reports. 

In the bright tail there were 165 observations. 
Light thunderstorms were reported in 29% of the observations 
with light rain being reported in 21% of the observations. 


Also, 19% of the observations were light rain showers and 


18% were moderate/heavy thunderstorms. Light precipitation 
was reported more often (72h) than moderate/heavy precipi- 
tation MWizsnye Convective type precipitation was also 


reported in 72% of the observations as opposed to stratiforn 
precipitation being reported in 28% of the observations. As 
expected, most of the moderate/heavy precipitation and thun- 
derstorm observations were in the bright end. 

In the region between the dim tail and the bright 
tail, there were 947 observations. Light rain Showers and 
light thunderstorms were dominate in this region with 30% 
and 29% of the reports respectiveiy. Light rain was 
reported in 22% of the observations. Light preci pi tation 
waS reported in 83% of the observations while moderate/heavy 
precipitation was reported in 1/% of the observations. 
Within this category, convective type precipitation was 
reported 74% of the time whiie stratiform precipitation was 
reported 26% of the time. 

An important consideration when evaluating the 


visual categories is the time of day. In this category, 57% 
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of the observations were within a three hour period fron 
Z200 to 90000 Giles This implies that most of the 
precipitation occurs late in the afternoon. Also, 37% of 
the 0.00 estimated albedo values (39) occurrec at 0000 GAT. 
Of the remaining 0.00 values, two occurred at 1300 GMT and 
one each at 2000 GMT, 2100 GMT, and 2300 GMT. This suggests 
that distortion caused some of the very low estimated albedo 


values. 


E. COMPARISON OF RESULTS TO PAUL (1983) 


ieee introduction 


This section will compare our results with the 
results by Paul (1983). Paul (1983) divided her data into 
seven precipitation categories: 

fe. -GOonvective: 

Pee COnt tnuous; 

fee light: 

YW. moderate/heavy; 

5. general; 

62 overcast ceiling; and 

7. overcast and broken ceiling; 
and two no-precipitation categories: 

i-ene-plCepitallVome Over CaSteycelling;: and 

Ze fo-~precipitation, broken and overcast ceiling. 
Each of the above nine categories were further divided into 
those with infrared satellite data and visible satellite 
data. 

As stated previously, this study used one precipita- 
tion category and one no-precipitation category, each 
divided into a visible and an infrared satellite group. Tn 
addition, only those surface observations which reported 


either broken or overcast ceilings were evaluated. 


a 


This chapter is divided into two Sections. One will 
compare results of the precipitation category and the other 


Will compare results of the rno-precipitation category. 


2s ) PECCIpItat ion 


— oa Set ee Oe 


This section will compare these results to those 
from Paul (1983) » for the precipitation with broken wer over 
cast ceilings (see Table VII). The mean values will be 


discussed first, followed by the standard deviations. 
a. MWeans 


The mean infrared values of this study are 
approximately 8.0 C colder than those fronmPaut ste (see 
Table VII). The differences can be attributed to the fact 
that nightime data are included in this study. In this area 
of the country many thunderstorms occur from late evening 
through the early morning hours. By including these storms 
(with their cold cloud top temperatures) in this study, the 
mean cloud top temperatures must decrease. 

Also, the means of the infrared data decreased 
slightly (became colder) in both studies as the matrix size 
decreased. This decrease can be attributed t9 the decrease 
in surrounding area when changing from a 10 x 10 matrix toa 
4 x 4 matrix. The smtailer matrix size allows the satellite 
to detect only the colder cloud» tops associated with the 
precipitation and not the warmer clouds from the fringe 
areas. 

There was very little difference in the visible 
mean values between the two studies (see Table VII). This 
1s not surprising because the time restrictions imposed by 
Paul (1983) contained most of the visible satellite data 
within the data set. The trend of the mean visual values is 
the same for both studies (see Table VII). Both show a 
Slight increase in albedo values when changing from a 10 x 


10 matrix toa 4 x 4 matrix. 
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b. Standard Deviations 


The infrared standard deviations are slightly 
Mreger in this Study than in Paul (1983) (see Table VII). 
The trend of increasing standard deviations with decreasing 
matrix size occurs in both studies. This can be related to 
the warmer fringe areas of the clouds being deleted by the 
smaller matrices. 

The standard deviations of the albedo values are 
danger in this study than in Paul (1983) (see Table VII). 
By restricting the time frame, Paul (1983) eliminated many 
possible distorted values that were included in the calcula- 
Prons for this study. Since the distorted values tend to 
accumulate at the extremes (0.00 and 1.20), including then 
in the calculations increases the standard deviations. 

The trend of the albedo standard deviation 
values is toward slightiy larger values as the matrix sizes 
decrease. The less fringe area covered by the 4 x 4 matrix 
diminishes the chances of the dimmer fringe area values 


being included in the calculations. 


3. No-precipitation 


ae Means 


The mean infrared values in this study are 
approximately 4.0 C colder than those from Paul (1983) (see 
Table VIII). The mean values of both studies remain almost 
constant as the matrix size is decreased. 

The visual means showed very little difference 
between the two studies (see Table VIII). As stated before, 
this is not surprising since most of the visual readings are 
within the time restrictions imposed by Paul (1983). In 
both studies, the mean albedo values tend to remain constant 


aoetne Matrix Size 1S decreased. 
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bz Standard Deviations 


The infrared standard deviations are about 1.5 C 
larger than Paul (1983) (see Table VIII). As in the 
precipitation category the standard deviations increase 
slightly as the matrix sizes decrease. 

The standard deviations of the visible values 
are approximately 0.06 larger than those in Paul (1983) (see 
Table VIII). This can be attributed to the addition of 
distorted mean values in this study which were exciuded by 
the time restrictions imposed by Paul (1983). The distorted 
values give additional weighting to each end of the range, 
causing the standard deviations to increase. The trend of 
the slightly increasing visual standard deviations with 


decreasing matrix sizes 1S again repeated here. 


rs SPADS THRESHOLD VERIFICATION 


This section will verify the no-precipitation/ 
precipitation threshold of a cloud and precipitation anal- 
ysis model of the United States Navy's interactive Satellite 
Data Frocessing and Display System (SPADS) using the 24 h 
data base. The objective is to present a comprehensive 
analysis of the infrared, visible, and infrared and visible 
thresholds combined. No attempt will be made to derive 
better thresholds, but only to verify the existing thresh- 
olas. After a prief explanation of the automated cloud and 
precipitation analysis program, there is a discussion on how 
thresholds work in an ideal Situation. Next is an explana- 
tion on how the verification was conducted, and then the 
results are presented. 

This program was designed to produce in real time (15 to 
30 minutes) analyses of important cloud and weather 
features. Te uses the Geostationary Observational 


Environmental Satellite (GOES) visual and infrared images to 
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mm@duce contoured digital displays o£ cloud types, cloud 
amounts, cloud top heights, cioud top temperatures, and 
precipitation intensity for an approximate 1024 x 1024 n mi 
sector of a geostationary image. 

The objective was to determine the success of the 
Peeecipitation/no~precipitation thresholds in separating 
Mieecifitation from no-precipitation for this data set. 
Ideally a threshold value will produce a clear and distirct 
separaticn between two events. Tf themeeevo distributions 
overlap (an area with eyual probability of either event 
Seeuring), the smaller the overlap the better. The 
@meeshOld with an overlap can still be useful as long as 
there is not too high a percentage of each event located 
within the overlap region. The more events iocated within 
the overlap region increases the chance of a misclassified 
event up to a point where there is ro longer any skill asso- 
ciated with this prediction technigue. 

The verification was conducted using the entire data set 
over the full 24 hours. The visual and infrared thresholds 
from the interactive cioud and pmecipitation analysis 
program were verified using the four categories descrived in 
Giapter II. The thresholds values used in the verification 
were; 

1. Estimated albedos greater than or equal to 0.55; and 
2. Cloud top temperatures less than or equal to 249.0 K; 
BOlmepLeCipitation, and: 
1. Estimated albedos less than 0.55; and 
2. Cloud top temperatures greater than 249.0 K; 
mone nO-precipitation. 

Using only the infrared threshold, the no-precipitation 
and precipitation histograms (see Fig. 2 and Fig. 5) were 
scanned for values which exceeded/were less than or eyjual to 
paosO K. The results (see Table IX) showed that 64% of the 
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precipitation observations had cloud top temperatures 
greater than or equal to 249.0 kK. Likewise, 87% of the 
no-precipitation observations had cloud top temperatures 
less than 249.0 K. 

Using only the visual thresholds, the no-»Drecipitation 
and precipitation histograms (see Fig. 3 and Fig. 6) were 
scanned for values which were greater than or equal to/less 
than 0.55. The results (see Tabie X) Showed that 63% of the 


precipitation observations had albedo values greater than 


Os 55. Likewise, 87% of the no-precipitation observations 
had albedo values less than 0.55. 
Next the above calculations were repeated, only this 


time the criterion for precipitation/no-precipitation was 
based on both the visual and infrared satellite data. If an 
observation did not contain both visual and infrared satel- 
lite data, it was not considered in our calcuiations. The 
resuits are shown in Table XI. They show that 82% of all 
precipitation observations in the data set had both albedo 
values greater than or equal to 0.55 and cloud top tempera- 
tures less than or equal to 249.0 K. Likewise, 93% of the 
ho~precipitation observations had albedos less than 0.55 and 
cloud top temperatures greater than 249.0 K. Using the 
visual/infrared thresholds together increases the chance of 
correctly classifying a precipitation event from 63% for 
visible and 64% for infrared to 82%. In addition, using the 
visual/infrared thresholds together increases the chance of 
correctly classifying a non-precipitation event from 87% for 


both visual and infrared to 93%. 


Ge SUMMARY 


This chapter consisted of analyses and discussions of 
our data set. It began with a discussion on choosing the 


most representative and reliable Matrix ize and the 


SZ 


distributions formed Ey those matrices. Nemiaetne G45 tia Du 
tions were tested for normality. Using those results, a 


discussion of the characteristics of the distributions were 


presented. From there, the results were compared to those 
Ouearned by Paul (1983). The final section was a veritfica- 
tion of the established precip? tatien/no-precipitation 


thresholds using the data set. 


23 


I¥V. SUMMARY AND CONCLUSIONS 


A. DATA DESCRIPTION SUMMARY 


A data set consisting of 70,623 surface observations, of 
whicn 29,342 had collocated satellite data, was studied. 
The satellite data consists of a 10 x 10 matrix centered on 
the surface observation station. The visual satellite data 
were converted to albedo values using the Muench and Keegan 
(1979) normalization scheme. This scheme also corrects for 
anisotropic scattering effects. The data set was divided 
into four categories of precipitation/no-precipitation and 
visual/infrared satellite data. It was further divided by 
uSing only those obServations which reported broken or over- 
cast ceilings. In addition, the full 24 hours of avaiiable 
data was analyzed, ccmpared to the 9 hours of daytime obser- 
vations used by Paul (1983). 


B. DATA ANALYSIS SUMMARY 


This chapter consisted of analyses and discussions of 
the data set. It was determined that the 4 x 4 matrix had 
the best combination of accuracy o£ location, and accuracy 
in best describiny the region based on a Singie observation. 
Next, adescription of how the histograms were forned and 
tested for normality was given. BecauSe one could not say 
with any acceptable degree of confidence that the distribu- 
tions were normally distributed, it was decided to analyze 
the non-normal parts and to gain a better understanding into 
the composition of the data set. The trends when moving 
from the cold/bright end to the warm/dim endS were as 
expected, but one could not accurately Separate light fron 


moderate/heavy precifitation uSing tne cloud top 
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Pwperdtures OF  albedos#when the results were compared to 
Paul (1983), there were many Similarities, especially when 
decreasing the size of the matrices from 10 x 10 down to 4 x 
4. 

The last area studied was the verification of estab- 
dished precipitation/no-precipitation thresholds used by the 
cloud and precipitation analysis scheme in SPADS with the 
data set. It was found that using both the visual and 
infrared satellite data together to determine precipitation/ 


no-precipitation produced the best results. 


C. CONCLUSIONS 


The conclusions of this study are: 

1. Based on the results of the Chi-Sguare test, the 
data set and all its subsets used in this study are not 
normally distributed; 

2. the histograms Or the precipitation/no- 
precipitation and infrared/visual categories each were 
distributed aS one would expect for precipitation and 
no-precipitation events respectively; 

See ALthnough the peaks of the pEeCcipltation/ne- 
precipitation and infrared/visual distributions were as 
Smected, the precipitation and no-precipitation distribu- 
tions display partial overlap; 

4. Investigation of the overlap areas indicates that 
Multiple and thick high cloud layers are responsible for 
erroneous estimates of precipitation within the 
no-precipitation distributions. Concerning the precipita- 
tion distributions, the scattered character of the summer- 
time precipitation MJlikely produced a significant number of 
misclassifications; 

5. Albedo estimates for the brightness normaiization 


scheme are not reliable in the low sun angle situations; and 
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6. Estimates of reliability £06 precipitation/no- 
precipitation threshclds were obtained from the “data seu 
Using both the infrared and visual precipitation/no- 
precipitation thresholds combined greatly improves’ the 
precipitation estimation when using each threshold sepa- 


rately. 


De. SUGGESTED FURTHER STUDY 


The recommendations for further study are: 

1. Further analyses of precipitation categories should 
be conducted by separating them into specific precipitation 
types (i.e. thunderstorms, rain showers, etc); 

2 eLo additional studies On precipitation/no- 
precipitation threshclds for cloud/precipitation estimation; 
and 

3. Gather and study additional coliocated data sets in 


other seasons and for other regions. 
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PREVIOUS RESEARCH BY PAUL (1983) 


A. BI-SPECTRAL AND INFRARED THRESHOLD 


Paijas (198ia, 1981b) developed a bi-sSpectrai cloud 
classification based on viSuai and infrared data from the 
polar orbiting TIROS-6 satellite (see Fig. 7). The data set 
consisted of a limited number of daily observations, chosen 
GBerpeetheir synoptic characteristics, in May and August 1979 
Over a region encompassing Norway, Sweden, Finland, and the 
Baltic Sea with weather charts providing the ground truth. 


Based upon the precipitation threshold results of Muench and 


Keegan (1979), Liljas chose a cloud top temperature 
mmesnold of —-i2Z2 C to -15 C to classify cumulonimbus and 
hnimbostratus clouds. Seeds Guam it teeth Sec loudaclassi fica— 


tion and the assumption that the highest and the densest 
clouds produce the maximum precipitation amount, Daniega.s 
Suggested a qualitative precipitation scale based on the sun 
of the visual and infrared satellite digital counts (see 
Wale XII). These sums represent the areas of the Liljas 
nimbostratus and cumuionimbus cloud types in his bi-spectral 
G€ioud classification (see Fig. 8). 

Mmeovejoy and Austin {1979) studied rain mapping of cloud 
areas based on GOES visual and infrared satellite data over 
Montreal, Canada, and the tPeopledi. “Atlantic (Global 
Atmospheric Research Program Atlantic Tropical Experiment, 
eee data) With radar providing the ground truth. The 
Montreal data set consisted of 17 obServations over three 
days during June 1977. Working with 4 x 4 resolution satel- 
lite image, Lovejoy and Austin plotted two dimensional 


Irequency grids for the radar-determined rain and no-rain 
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points on a 25 x 25 array "(see fig. 9 andra the 
visual data were normalized by selecting the "brightest" and 
"dimmest" values in each image and iinearly interpolating 
the radiances between 0 and 1. 

Lovejoy and Austin (1979) state, with reference to the 


cumulus rain data distribution of Fige 9 that iesdiseem 


butl OnmewaS a) tO ued good approximation a two-dimensional 
Guassian." They do not provide or describe the statistics 
to support this assertion. The no-rain cumulus cases (Fig. 


10) were described as a bimodal distribution with one peak 
near the low visual and low infrared vaiues and the other 
peak near the rain peak but shifted slightly toward lower 
values, In most cases, the separation of the cumulus rain 
and no-rain caseS waS StatiStically Significant with the 
probability ranging from 10% to 50% that the rain and 
no-rain samples came from the same popuiation. 

The Lovejoy and Austin (1979) two dimensional frequency 
plots for non-cumulus storms were limited to one case. The 
Significant differences between the cumulus and non-cumuius 
data sets were the non-cumulus no-rain plot lost its bimodal 
character, relative to the cumuluS no-raia plot, and 
appeared as a broad two dimensional Gaussian distribution. 
The non-cumulusS rain plot points fell within the no-rain 
distraputaou- but were shifted slightly higher in the 
visual. The separation of the non-cumulus rain and no-rain 
cases was not statistically significant, with greater than a 
50% probability of the rain and no~-rain samples coming fron 
the same population. 

Lovejoy and Austin (1979) attempted to further classify 


the cumulus rain and no-rain cases into no-rain, light rain, 


and heavy rain. Rainfall cates greater than 2 mm-h-1, as 
determined by radar, were defined as heavy frain. As 
expected, the mean of the heavy rain cases was’ shifted 


Siightly toward higher visual and infrared values than the 
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mean of the light rain cases. However, the shift was so 
small that there was at least an 80% orobability of the 
light rain and the heavy rain caseS coming from the same 
population. Lovejoy and Austin (1979) concluded that 
"little, if any rainfall-rate information 1S contained in a 
Single (visual and infrared) satellite image. 

Lovejoy and Austin (1979) tested a spectral threshold 
technique for rain area mapping. Each satellite image of 
400 x 400 km was divided into one hundred 40 x 40 km boxes. 
The 100 subareas were each checked with radar to determine 
the total number of rain areas. An egual total number of 
Satellite subareas were classified as fain areas. The 
Satellite subareas with the highest visual and highest 
(cold) infrared values were classified aS rain areas, until 
the total number of satellite rain areas equaled the total 
number of radar determined rain areas. Phy SmeS pect Fa! 
threshold technique was applied to three days accumulation 
Geedata and 1S Shown in Tables XIII and XIV. When comparec 
to the success of the two dimensional freyuency plot method, 
the visibie and infrared thresholds averaged 45% and 58% 
worse, respectively. iiemeacctracy eo pthe visual threshold 
is limited by the extent of low, thick clouds and the 
infrared threshold 1s limited by tne extent of the cirrus 
clouds in the satellite image. Lovejoy and Austin (1979) 
concluded that "the errors involved in uSing a 'best thresh- 
old' are very large indeed." 

Del Beato (1981) studied correlations between cloud top 
temperatures (based on NOAA-5 satellite data) and rainfall 
teorals for 30 and 60 minute intervals Over eastern 
Ast ralla. The satellite data had a 60 sy km maximum reso- 
lution at subsatellite point and cloud top temperatures were 
area averaged for resolution of 200 sq km. The 21 data sets 
were first classified according to synoptic Situation ina 


rough attempt to group the data by cloud type, droplet 
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spectra, and alr Mass (tranceeon,- The initial results 
suggested that the cloud top temperature determined an upfer 
limit on rainfall amount, with the maximum increasing as the 
cloud top temperature decreases. A linear correlation anal- 
ysis to determine a quantitative relatiorship between rain- 
Tall amount and cloud top temperature gave indefinite 
results. 

Further study of surface and radiosonde observations 
indicated that classification by proportion of cunulitomn 
Cloud reports to allcloud reports and subcioud layer 
humidity might be more appropriate (Del Beato, 1981). This 
classification resuited ina correlation coefficient of 
0.90, excluding cases with cumuliform portions less than 50%. 
and dew point depressions of greater than 6 C. Fin diijyeee 
composite frequency distribution was calculated based on 
three cases, all southwesterly stream Situations described 


as "post-frontal cellular convection cases in cyclonically 


curved flow." The fitted equation was: 
£ = 0.057 —- O200GeT 1 —F Ones. (A.1) 
where f is the rainfall frequency, R is the 30 minute rain 


total (mm), and CTT is the cloud top temperature in degrees 
Celsius. The equation was fitted to 41 independent f 
values. This equation is associated with a correlation of 
0.79 at the 99% confidence level. Equation A.1 indicates no 
Cain from clouds warmer than +13 5 and a maxinum 30 minute 
rainfall of 2.5 mm for a cloud top™tenperacure sot —Z0ee 

Ih summary, Del Beato (1981) found that cioud top temp- 
eratures and 30 and 60 minute rainfall totals jandicated 
statistically significant relationships for cloud systems 
with a high portion of cumulus clouds and high subcloud 
bum d ea Additionally, as cloud top temperatures decrease 


to at least -35 C, rainfall totals increase. 
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Wylie (1982) attempted to correlate rainfall occurrence 
with radiosonde soundings, hourly Service-A observations, 
and visual and infrared satellite data. dis data sample was 
restricted to "large-scale cloud cover" areas with wide- 
Spread precipitation (rain gauge reports varied less than 
20%) for the Great Plains States region for the period 27 
February 1981 througo 4 January 1982. From 13 parameters 
derived from three data sources (see Table XV) the best 


linear regression equation for estimating rainfall rate was: 
Pieerain = 1.0242 + 0.380Pw - 0.0304Qc - 0.0047Ct (A.2) 


where Pw is the vertically integrated precipitable water 
vapor (in), Qc is the moisture convergence (g/kg/day), and 
Ct is the cloud top temperature (Kelvin). Equa tion A. 2ehas 
a linear correlation coefficient of 0.60. Linear regression 
eguations were also determined for the three parameters 
alone and for a combination of Pw and Qc to be used when not 
all three data types were available. The cloud temperature 


regression eguation was: 
Paeraaem (inj = 2.10 - 0.008Ct (A%'3) 


The correlation coefficient was -0.35. Wylie (1982) stated 
that toe synoptic scale data base measurements were best 
Suited for estimating broad changes in rainfall rates asso- 
Cliated with air masses and not suited for eStimating rain- 


fall rates associated with small scale dynamic processes. 


Be. LIFE HISTORY 


The life history methods are empirically derived precip- 
itation estimation schemes based upon two assumptions, 


PLES t, that significant rainfall cones frrom convective 


4] 


clouds, and second, that convective clouds can be identified 
and measured in satellite images. These methods involve 
Manual analyses of convective cloud areaS in a Sequence of 
visual, infrared, or both visual and infrared satellite 
images. Threshold values and study condition parameters 
associated with published life history studies are summa- 
rizedwen Table XVI. 

The Scofield/Oliver (Scofield, 1981) analysis follows a 
cecision tree procedure to estimate half-hourly rainfali for 
deep convective systems within tropical air masses. Using 
enhanced infrared and high resolution visual satellite 
images, the technique involves first identifying the active 
convective portion of the cloud, or cluster, from two 
consecutive satellite images. Once the active portion is 
identified, the halilf-hourly rainfall estimation is computed 
baseaq on such factors as cloud top temperature, cloud 
growth, and departure of precipitable water from a summer- 
time normal. 

The Griffith/Woodley (Griffith et al, 1978) technique is 
designed to estimate rainfall in the tropics, over large 
Space and time scales, using geosynchronous visual or 
infrared satellite imagery. This time-dependent technique 
waS empirically derived as a relationship between cloud 
area, echo area, and rain rate for two areas in south 
Florida, with» raingauge=radar™providimg the ground trie 
and was then tested in other tropical areas. This scheme 
was subsequently tested further in extratropical areas 
(Griffith et al, 1980), with modifications to the painitage 
amount predicted. 

The determination of a cloud area-rainfall relationship 
first required the specification of both a visual and an 
infrared threshold to define the cloud area. The visual 
brightness threshold, normalized for radiation geometry, was 


80 counts for the third Application Technology Satellite 
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ees -3) and the ianttared threshold was 253 K (-20 C). The 
thresholds were based on a comparison of the clouds with a 
given maximum digital count and the radar echoes associated 
with these clouds. 

The empirical cloud area-rainfall relationship was 
derived aS a two step process. First, a relationship 
between tne cloud area and the radar echo area, normalized 
for the maximum area achieved by the cloud or the cluster, 
was established for the visible and infrared satellite data. 
Second, the relationship between the echo area and the rain 


volume was determined and was of the forn: 


Rv = jAe (A.4) 


where Rv is the rain voiume per hour (m3 ~-h-1), I iS rain in 
iets Of (m3 ~km-2 —h~1), and Ae is tne echo area (kn2) 
defined py the Imm -h-1 rain rate. Thus, given a tine 
sequence of convective clouds (or cluster areas) measured 
from visible or infrared satellite images, a volumetric rain 
rate can be estimated. 

Seoutceet al, (1979) modified the Griffith/woodley tech- 
higue (Griffith et al, 1978) to estimate volumetric rain 
rate directly from a cumulonimbus cloud area and area change 


meeGording to the equation: 


hee—ea (0)A + a(1)dA/dt (Ave) 


where R is the volumetric rainfall of the cloud (m3 -s-1), A 
is the cloud area (m2), dA/dt is the change of cloud area 
over time (m2 ~s~1), anda(Q) and a(t) are constants with 
dimensions (m ~-S-1) and (m) respectively. The two constants 
were calculated by a least square fit of cloud area-rain 
rate pairs based on visible and infrared geosynchronous 


satellite data and 5.3 cm Ship radar rain data coiiected 
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during GATE. The cloud area and its change are defined by 
the threshold value. The visible threshold for cloud area 
Calculations was 60 digital counts on the ATS-3 £4(corre- 
sponding to an albedo of 0.45 with the sun overhead), or 
equivalently 172 digital counts on the Geosynchronous 
Meteorological Satellite (SMS 1). The infrared threshold 
was, 160 digital counts (—26.¢) The standard error between 
the estimated rainfall and the mean radar rainfall was 62% 
and 76% for the visual and infrared eguations respectively. 

Wide (1979) attempted to use the tropical convective 
rainfall techniques of Griffith et al, (1978) and Stout et 
al, (1979) for estimating precipitation in Montreal, Canada. 
Using visual satellite data, corrected for the changing sun 
angle (Mosher, 1975), infrared satellite data, and 10.0 cm 
radar measured rainfall rates, Wylie studied six days of 
precipitation, three days each in June and September 1977. 
Wylie concluded that because of air mass differences between 
Montreal and the tropics, the Griffith and Stout estimation 
techniques did poorly in Montreai, Canada. The singlemost 
important limitation with these two schemes was tke diffi- 
culty in measuring cumulonimbus cloud area when the "anvils 
were often merged into large cloud masses and the extensive 
Stratus cloud cover often obscured the pictures." Wylie 
also noted that the Griffith et al, (1978) threshold of -26 
C had to be changed to -16 C for the summertime Montreal, 
Canada area. With the warmer cloud top temperatures the 
cloud area were a larger, more appropriate size for 
tracking. 

Wylie (1979) then attempted to combine sounding data 
input into a one-dimensional model (Simpson and Wiggert, 
1969) and satellite cloud cover measurements to estimate 
rainfall for Montreal. With the GATE measurements for rain 
rates associated with satellite- derived cloud areas and the 


model output, rainfall rates were estimated by multiplying 
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the two values. The most accurate eStimations were for the 
Cumulus clouds in the warM air masses occurring in June, the 
cases the model was designed to handle. Wylie concluded 
that in crder to estimate the rainfall in all geographical 
areaS and seasons amore sophisticated model would be 
needed. 

Negri and Adler (1981) did one case study of 15 thunder- 
storms in the Oklahoma, Arkansas, and 4issouri area on 24 
Meee! 1975. miey used radar data foreground truth and had 
special 5 minute GOES-E satellite passes over the area of 
interest. They were able to determine that precipitation 
began falling, aS indicated by radar data, for Cloud top 
Pemperatures ranging from 229 K to 260 K (-44 C to -13 C). 


The mean cloud top temperature value was 247 kK (-26 C). 
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APPENDIX B 
HUENCH AND KEEGAN (1979) NORMALIZATION SCHEME 


The Muench and Keegan (1979) normalization relates the 
normalized reflectivity (x£,) to the varying solaz@anglesama 
Maximum digital counts through the reflectance term (fF), and 
the anisotropic scattering through the X tern. Table XVII 
defines the symbols, Table XVIII lists the geometric iden- 
tity equations, and Table XIX lists the normalization equa- 
tionsce Fig. 11 gives an example of the normalization 


applied to the stated location. 
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APPENDIX D 


PRB LES 
Table i. Information Available for Each Observation 

le Station ID 
Zs Latitude 
Bie Longitude 
4. Altimeter 
Sr. Temperature 
Ge Dew Point 
7. Wind Direction 
8. Wind Speed 
or Sea Level Pressure 
10. Cloud Group | 
lt. WVisibilinge, 
12. Precipitation Amount 
13. Wind Gusts 
14. Current Weather 
15. Maximum Tenperature 
16. Mininum Temperature 
17. Snow 
18. Observation Date/Time 
19. Visible/Intrared Satellite Data 
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Table III. No-~Precipitation, Infrared Anageais 


Mean; 275.71 St Dev: ‘A927 
COLD TAIL: 1461 Obs ( < 240.0 K ) 
# Obs High Middle Low 
Cvercast - 0.50 S34 0.10 O05 
Broken 0.50 0.41 Ois0e7 O03 
Total 1.00 On76 0.16 as 


Table IV. No-Precipitation, Visible Analysis 
Mean: 2129 St Dev: ZOnriZ 

BRIGAT TAIL: 999 Obs (2>0e70 ) 

# Obs High Middle 

Overcast 0.49 O21 0.14 

Broken 0.51 0.41 0.07 

Total 100 0.59 Oe ee 
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Table VI. 


Precipitation, 


Mean: 242.99 
Cold 
End 
TRH O15 
TRW- Oe 
RW sro 
RW- 0.20 
R 02,02 
R- 0.14 
TR 0.00 
TR- 0.05 
# Obs 959 
PEcei pi tation, 

Beams. 61432 
Bie 
End 
TRW 0.18 
TRW= O29 
RW 0.06 
RW- 0.19 
R 0202 
R- 0.21 
TR 0.01 
TR- 02°03 
# Obs 163 
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Table VII. Comparison of Means and Standard Deviations for 
Precipitation with Broken/Overcast Ceilings 
with Paul { 1l2ca 72 (Paul's values in 
parenthesis) 


Precipitation - Infrared 


cme 3 4x 4 
Mean 244.0 Kee (252-6 243.0 KR” (250.6 5% 
St Dev 21. Tet. 20. Sa 21.7 Kk ( 25 
Precipitatiou - Visible 
10 x 10 4x 4 
Mean 0. 567m 0. oo) 020155) (Ceo 
St Dev Os-2 1 Ze lee) 0.285 (02225) 


Table VIII. Comparison of Means and Standard Deviations for 
No-precipitation with Broken/Overcast Ceilings 


with Paul (1983). (Paul's values in 
parenthesis) 
No-Precipitation - Infrared 
10 x 10 ex 
Mean 2175-0 Ka 2 7 Seem 213560 Ko (2752 ee 
St Dev 18.4- Kk ( 16260) 19.4% ( 172858 
No-Precipitation - Visible 
10 x 19 4x4 
Mean On ZO Sm Oz) O22 73 (Om27a) 
St Dev O..250 9 O22) Oe 26305 ( 007) 
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Table IX. VeomrrELeatilonmormesPADS.~Lafrared. Thresholds 


Observations: Yes No om au 
satellite fe. 7 — - Obs 
Estimate: 7 

Yes 1195 (0.64) 661 (0.36) 18 56 

No 2503 (0.13) 167851 (0.37) 19,354 

Total Obs 3696 ponies 2A lO 
Taole X. Verification of SPADS Visible Thresholds 

.. Observations: Yes No Total 
Satellite ~ Obs 
Estimate: 

Yes 309 810263) Be Ot Ole30 ) 1437 
No 1962 (0.13) UZ, 78 (0-37) 14,740 
Total Obs 2871 13306 16,177 
Hable XI. Verification of SPADS Visible and Infrared 
Thresholds Combined 

Observations: Yes No Total 
Sieeilite . a - Obs 
Estimate 

Yes 843 (0.82) so) (ogee) 10 23 
No 825) (0207) 107 ao (0.93) 11,144 
Total Obs 1668 10,499 V2 lo 
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Table X£27. 


Table XIII. 


Threshold Values eee Precipitation 
Intensity Levels as Applied in Fig. 8 
(Liijas, 1981a) 
The Sum of Digital Levels 

Chl..+-Gheu 291 - 310 Ligh taka? ni 

Sle) 330 

2) — on 

35a = 370 

si = 390 

39 3+ 


Very Strong warn 


Statistical Comparison of Rain Area Mapping 


Technigues (Lovejoy and Austin, 1979) 
Opt. 2-D Boundary IR Optimum Threshold Visible Optimum Threshold 
2 ee a Area Rain 
Day (Ra/R) x 100° IR(K) (Ra/R) x 100 (Scale: 0-1) (Re/R) x 100 Coverage (7) 

GATE 
242, 243, 246 65 < 232 53 >0.68 64 15.8 
247, 248, 251 
252, 261 
Montreal 
152 56 < 232 20 > 0.88 31 9.7 
180 56 < 247 55 >0.80 48 24.0 
153 53 < 254 52 >0.88 49 15.9 


*Referred to in text as “percentage of correct satellite rain.” 


Taple XIV. 


Technique 


2-D Pattern 
Matching 

2-D Pattern 
Matching 

Optimum IR 
Threshold 

Optimum Visible 
Threshold 


Statistical Comparison of the Accuracy of 
Areas (Lovejoy and Austin, 1979) 


Total No. 
of Points 


47706 


40361 
33738 
22558 


Rain 


Eams 


0°22 
0.18 
0.71 


0.58 


Number of Images Error 

Region or Sequences Bias Factor 
Montreal 17 la 1.26 
Montreal 3 1.08 1.19 
Montreal 3 1.38 1.74 
Montreal 3 1.54 le? 
GATE 8 te21 1.41 


2-D Pattern 


O75 


Matching 


Zz 


Table XV. General Coefficients for the Determination of 


Precipitation Based 


(Wylie, 1982) 


GORR ELATION 
PARAMETS 8 PRECIP. REPOR? 
epee ci ty integrated 0.48 
precipitacle water vapor 
@ioud top brightness “0.44 
Cloud top height =()2.40 
Moisture convergence 02338 
Cloud top temperature ners © 
Bubble model predicted cond. Oe 
500 mb vorticity advection -0.21 
Parcel lifted index =). 20 
700 mk temperature advection 042.20 
Sfc temperature advection 0.19 
850 mb temperature advection Ocean: 
Wind convergence (sfc) 0.09 
Vertical wind shear 0203 


wor acant ccrrelation at 


the 99% lsavel. 
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184 


190 


184 


199 


115 


Tis 


200 


173 


156 


189 


167 


156 


on the Three Data Types 


try 


ao 


44x 
Sror* 
on™ 


eure 


8 x 
8* 


7 * 


Table XVI. 


STUDY 


Griffith, 
et al., 
(1978) 


Stout, 
et al., 
(1979) 


Wylie 
(1979) 


Negri and 
Adler 
(1981) 


Summary of Life BRistory Threshold Values 


LOCATION 


Florida, 
Venezuela, 
Honduras, 

and hurricanes 
impacting East 
Coast United 
States 


tropical North 
Atlantic 


Montreal 


Oklahoma, 
Arkansas, 
Missouri 


* ATS=-3 satellite 


CASES 


34 days 


57 obser-— 
vations 


6 days 


l day 


(15 thunder- 


storms) 


TIME OF 
YEAR 


summers 
1969-1976 


September 
1974 


June 1977 
September 
1977 


April 24, 
1975 
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THRESHOLD VALUE 


INFRARED 


S056 


~26°C 


=T626 


27°C 


VISUAL 


x 
80 counts 


0.45 albedo 


(sun over- 
head) 


moapie XVII. 


Symbol 
C 
C 


Pest sOL 2a 1 bOLS 


Description 
GOES video count number (0-63) 


GOES video count number for perfect 
diffuse reflector and overhead sun 


Greenwich meridian time 


Distance of earth to sun 

Mean distance of earth to sun 
Julian date 

Hour angle 

Cloud reflectivity 


Arc-length observer to subsatellite 
jeje pele 


Declination of the sun 

zenith angle of the sun 

Longitude 

Longitude of subsatellite point 
Anisotropic scattering coefficient 
Latitude 

Azimuth of the sun 


Azimuth of the satellite 
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Units 
dimensionless 
dimensionless 
hours-minutes- 
seconds 
kn 
km 
dimensionless 
radians 
dimensionless 


radians 


radians 
radians 
radians 
radians 
dimensionless 
radians 
radians 


radians 


Table XVIII. Basic Geometric Satellite-Earth Relationships * 


Declination: 

6 = 0.408 sin [(d-81) * 27/365] 
Solar distance ratio: 

R/R, = 1 - 0.167 cos [(d-14) * 217/365] 
Hour angle: 

h = A + mw = G(hours) x 1/12 
Arc-length: 

cosy = cos(h. - A) cos¢ 
Satellite azimuth: 

sin ($5 - *) = sin (Ag - A)/siny 
Solar azimuth angle: 

cos¢ = siné siné + cos cosé cosh 
Solar azimuth: 


Sing, = cosé sinh/sine 


* 
Angles in radians 


Table x Muench and Keegan (1979) Normalization 
EgGuations * 


Ga.2 
r = (=) * secl 
Cy 
4¢ = |o5 - $4 
2 
C, acosme (€(F —- 50) "= 178) 
Co wm O027 COS ((C = 22.5) *s4p ee cose ) 


Ce - Sone ((4d - 70) * 1.3) 


Dad 
We 


1.0 + 0.05 * (1 + cos(2*z)) + 0.20 * (Cy 7 Co) . C. 


~- 


ry [1.09 - 2*(1.09 - F * y *(R/R,)*)/(1 + cost/? Ei] 


*« 
Angles in degrees 
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